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ABS TRAC T 

The feasibility of using active controls to delay 
the onset of whirl-flutter on a joined-wing tilt- 
rotor aircraft was investigated. The CAMRAD/JA 
code was used to obtain a set of linear 
differential equations, which describe the 
motion of the joined-wing tilt-rotor aircraft. 
The hub motions due to wing/body motion is a 
standard input to CAMRAD/JA and were obtained 
from a structural dynamics model of a 
representative joined-wing tilt-rotor aircraft. 
The CAMRAD/JA output, consisting of the open- 
loop system matrices, and the airframe free- 
vibration motion were input to a separate 
program, which performed the closed-loop, 
active control calculations. An eigenvalue 
analysis was performed to determine the flutter 
stability of both open- and closed-loop systems. 
Sensor models, based upon the feedback of pure 
state variables and based upon hub-mounted 
sensors, providing physically measurable 
accelerations, were evaluated. It was shown 
that the onset of tilt-rotor whirl-flutter could 
be delayed from 240 to above 270 knots by 
feeding back vertical and span-wise 
accelerations, measured at the rotor hub, to the 
longitudinal cyclic pitch. Time response 
calculations at a 270-knot cruise condition 
showed an active cyclic pitch control level of 
0.009 deg, which equates to a very acceptable 9- 
pound active-control force applied at the rotor 
hub. 
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NOMENCLATURE 

Aj (i=0,1,2) system matrix 
A first-order system matrix 

Bo control matrix 

B first-order control matrix 
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Indices: 

a 

anti-symmetric mode 

ac 

active control 

eg 

center of gravity 

F 

airframe body axis system 

FT 

trim Euler angles 

P 

pilot control, or airframe location P 

h 

hub 

h,c 

hub motion , chord-wise direction 

h,s 

hub motion , span-wise direction 

h,v 

hub motion, vertical direction 

s 

symmetric mode 
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INTRODUCTION 

The goal of the tilt-rotor concept is to achieve 
the cruise speed of a fixed-wing aircraft, while 
retaining the hover capability of a helicopter. 
Current tilt-rotor aircraft (XV-15 and V-22) 
employ very thick wing airfoils (23% thick) to 
obtain adequate wing stiffness and strength to 
handle the loads imposed by vertical-jump 
takeoffs and high speed whirl-flutter stability. 
These thick airfoils result in wing 
compressibility effects limiting the high speed 
potential of the tilt-rotor aircraft, e.g. to 
M-0.575 for a XV-15 size aircraft. In hover 
flight the net thrust of the tilt-rotor aircraft is 
reduced by the vertical drag (download) of the 
wing 1 . Wolkovitch, et. al. 2 , performed an 
analytical study to evaluate the application of 
the joined-wing concept 3 to tilt-rotor aircraft. 
Figure 1 shows the baseline cantilever-wing 
tilt-rotor configuration, representing a XV-15 
size tilt-rotor aircraft. Figure 2 shows a 
typical joined-wing tilt-rotor configuration 
studied in Ref. 2. The joined wing concept can 
reduce the projected wing area in the rotor 
downwash in hover, thus potentially reducing the 
hover download. The joined wing concept also 
allows for thinner airfoil sections (12% thick), 
resulting in an increase in the drag-divergence 
limiting Mach number from 0.575 for the 
cantilever-baseline model to 0.69 for the joined- 
wing model. 

Reference 2 evaluated various joined-wing tilt- 
rotor configurations in terms of potential 
performance improvements, airframe aeroelastic 
characteristics, and aircraft whirl flutter 
speeds. Design parameters that were varied 
included front wing sweep angle, front wing 
airfoil thickness, and nacelle center of gravity 
location. 

Reference 2 showed a potential speed 
increase of 100 knots by delaying the 
compressibility drag effects to higher Mach 
numbers. However, this performance advantage 
was negated by a lowering of the joined-wing 
whirl-flutter speed from 330 to 240 knots. 
Representative results for Ref. 2 are shown in 
Fig. 3. 

Studies have been performed in the use of active 
controls to reduce a tilt-rotor's response to a 
gust 4 ' 7 and to reduce tilt-rotor blade loads 
during maneuvers 8 . Other studies have been 
performed to improve the hover aeromechanical 
stability of helicopters 9 . Clearly an area of 
interest for tilt-rotors is flutter control with 
the objective to allow the tilt-rotor aircraft to 


fly at higher speeds or to reduce its structural 
weight. 

Nasu 10 studied flutter control on a simplified 
tilt-rotor model, consisting of a semi-span 
straight wing, a pylon attached to the wing tip, 
and a three-bladed hingeless rotor. Nasu 

considered symmetric wing elastic motion (wing 
beamwise, chordwise, and torsional bending) and 
rotor blade chordwise and flapwise bending. His 
feedback control was performed in the state- 
space domain. The author studied the use of 
active controls to delay the aeroelastic 
instability (whirl-flutter) on a tilt-rotor model, 
representative of a XV-15 size aircraft 11 . 

The results of a whirl-flutter alleviation study 
for a joined-wing tilt-rotor aircraft are 

presented here. Table 1 provides a comparison 
of the tilt-rotor models studied in Refs. 10 and 
1 1 with the present investigation. This study 
included rigid-body motion, both symmetric and 
anti-symmetric elastic body degrees of freedom, 
rotor blade bending and torsion modes, the hub 
gimbal and the rotor speed degrees of freedom. 
In addition, feedback control in the state-space 
domain as well as feedback control using 
realistic sensor models (the output of hub- 
mounted accelerometers) was investigated. 

INVESTIGATION APPROACH 

The approach followed in this joined-wing tilt- 
rotor whirl flutter alleviation study is described 
in this section and is schematically shown in 
Fig. 4. 

A number of XV-15 size joined-wing tilt-rotor 
configurations were studied in Ref. 2, which 
used the structural analysis code MSC-PAL to 
calculate the airframe mode shapes, modal 
frequency, and modal mass. Only the airframe 
mode shape at the hub is reported in Ref. 2. 

A typical joined-wing tilt-rotor aircraft 
configuration (model 166CL in Ref. 2) was 
selected to study the potential of using active 
controls to delay the occurrence of whirl flutter 
on a joined-wing tilt-rotor. The airframe 
structural dynamic characteristics, consisting 
of the modal frequency, modal mass, and the hub 
motion due to elastic airframe deformation 
were obtained from Ref. 2 and are used as an 
input to the rotorcraft analysis code 
CAMRAD/JA, which provided the mathematical 
plant model of the tilt-rotor aircraft. The 
CAMRAD/JA analysis was also used in Ref. 2 to 
determine the whirl flutter speed for the various 
joined-wing tilt-rotor aircraft configurations. 
The CAMRAD/JA (Comprehensive Analytical 


Model of Rotorcraft Aerodynamics and Dynamics, 
Johnson Aeronautics) is described in detail in 
Refs. 12-15 and a summary description is 
provided in Appendix A. 

CAMRAD/JA is used to develop the equations of 
motion for the airframe motion, the equations of 
motion for the rotor, and the expressions for the 
rotor hub load reactions. The coupled airframe- 
rotor equations of motion are obtained by 
substituting the hub motion into the equations 
for the rotor motion and hub loads, and then 
substituting the hub reactions into the airframe 
equations of motion. The aircraft motion 
consists of the six rigid body degrees of freedom 
and the airframe elastic free-vibration modes. 
CAMRAD/JA accounts for the rigid body motion in 
a quasi-static manner. The CAMRAD/JA analysis 
specifically uses the MSC-PAL data in modeling 
the airframe dynamics at the rotor hub location 
due to the elastic airframe deformation. 

For purposes of an aeroelastic analysis, 
CAMRAD/JA linearizes the aerodynamic and 
inertial forces around the aircraft trim solution 
to derive a mathematical plant model, consisting 
of a set of linear differential equations 
describing the perturbed motion of the aircraft 
about the trim condition. 

The CAMRAD/JA plant model was used as an 
input to a separate program, which allowed for 
the evaluation of various feedback control 
system schemes for their effectiveness in 
delaying the occurrence of whirl-flutter. 
Additional input to the closed-loop analysis 
program consisted of the motion sensor model 
and the active control feedback specifications. 
The analysis code uses an eigenvalue analysis to 
determine the stability of the closed-loop 
system and allows for the calculation of the 
system time-history response to various types 
of control inputs, such as an impulse or step 
control input. 

A separate utility was used to obtain a 
mathematical formulation describing the output 
of the airframe-mounted motion sensors. Input 
to this utility consisted of the geometric 
definition of the sensor location on the airframe, 
the airframe mode shape data for the sensor 
location, and CAMRAD/JA parameters describing 
the aircraft trim conditions (i.e., the aircraft 
velocity and orientation, and its flight path). 

Sensor models, based on pure state variable 
domain feedback as well as based on realistic 
sensors, measuring the actual "physical" motion 
of the airframe at the hub were evaluated. An 
eigenvalue analysis was used to determine the 


stability of the aircraft with the closed-lobp 
feedback system. Time response calculations 
were performed to determine the 
magnitude/level of the required active control 
inputs. The effect of sensor noise on the closed- 
loop aircraft stability was also investigated 
using random noise of a specified magnitude to 
simulate measurement noise and/or system 
modelling uncertainty. 

THEORETICAL FORMULATION 

PLANT MODEL DEVELOPMENT 

For the purpose of performing an aeroelastic 
stability analysis, the rotorcraft motion can be 
described by means of a set of linear differential 
equations of the form 12 : 

A 2 X + Aix + Aox = Bov (1) 

where x is the vector of degrees of freedom, v is 
the vector of controls, Aj (i=0,1,2) are the 
system matrices, and Bo is the control matrix. 
Equation (1) represents the rotorcraft plant 
model, i.e., the tilt-rotor aircraft without 
feedback control, referred to hereafter as the 
basic aircraft. The Aj (i=0,1,2) and Bo matrices 
are obtained from the CAMRAD/JA code and are a 
function of the cruise flight conditions. For the 
most general case in which Eqn. (1) describes 
the motion of a rotorcraft, the matrices Aj 
(i=0,1,2) and Bo have periodic coefficients and 
must be analyzed for stability by the methods of 
Floquet-Liapunov theory. For the tilt-rotor in 
cruise flight, the rotor operates in a mostly 
axial-flow environment and the Aj (i=0,1,2) and 
Bo matrices in Eqn. (1) are constant coefficient 
matrices. Consequently, the techniques for the 
analysis of time-invariant systems can be 
used 16 . 

The CAMRAD/JA basic aircraft plant model 
matrices (Eqn. (1)) are input to the closed-loop 
analysis code which models the airframe motion 
sensor outputs and feedback control system as 
described by Eqns. (2)-(8) below. 

The control vector v in Eqn. (1) consists of the 
sum of the pilot control inputs, v p , and the active 
control inputs, v a c: 

v(t) = v p (t) + vac(t) (2) 

The motion sensor output is described by the 
observer or sensor vector, y, which is defined 
as: 
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y * C2X + Cix + Cox + Dov (3) 

where Cj (i»0,1,2) are the system sensor 
matrices and Do is the control sensor matrix. 

The active control vector, vac. is defined by: 

Vac = Fy (4) 

where F is the output feedback gain matrix. No 
attempt is made here to model the dynamics of 
the swashplate actuators. 

Combining Eqns. (3) and (4) and substituting 
(Vp+Vac) for v in Eqn. (1) leads to the following 
equation for the closed-loop system: 

(A2 -BoFC 2 )x + (Ai-BoFC-|)x + 

(Ao-BoFCo)x = (Bo+BoFDo)v (5) 

which can be rewritten as: 

A|x + A*ix + A ox = Bjv (6) 

Equation (6) is transformed to a first-order 
state variable form by the substitution: 



yielding: 

x s = Ax s + Bv (8) 

An eigenvalue analysis can be performed to 
determine the stability of the rotorcraft system 
represented in the state-space form by Eqn. (8). 
The natural frequency, co n , and the damping ratio, 
C, (fraction of critical damping), for each mode is 
obtained from the associated eigenvalue, X, as 
(o n -|X|n/2n and £=-Real(X)/|X|, respectively, 
where ft is the rotor rpm in rad/sec. 

SENSOR MODEL 

The sensor model used in the feedback control 
system is defined by Eqn. (3). By appropriate 
selection of the elements of the Cj (i=0,1,2) by 
setting the elements to ”1* or *0”, a feedback 
scheme based upon pure state variables can be 
obtained. Alternatively, the output of a sensor, 
such as an accelerometer mounted on the 
airframe can be modeled by the proper 
specification of the elements of the Cj matrices. 
The coefficients in the j ,h row of the Cj matrices 
represent the contribution of each aircraft 


degree of freedom to the sensor output signal 
represented by the corresponding j th element of 
the sensor vector y. 

The contribution of k th aircraft degree of 
freedom to the output of a sensor located at a 
point P can be the due to the local airframe 
displacement, x Pi k, (through the corresponding 
Co-matrix element), the velocity at point P, k p k 
(Ci -matrix element), or its acceleration X Pi k 
(C2-matrix element). 

It is assumed that single-direction motion 
sensors are employed in the sensor model. 
Therefore each element of the sensor vector y 
represents the motion of point P in only one 
direction. Only airframe-mounted sensors are 
considered. The rotor degrees of freedom do not 
contribute to the sensor output. 

To determine the element values of the Cj- 
matrices (i=0,1,2) it is necessary to determine 
the aircraft body motion at point P where the 
sensor is located. The derivation for this body 
motion at point P follows the derivation for the 
hub motion as presented in Ref. 12 and is 
described below. 

The body axis system is defined as pictured in 
Fig. 5. From the pilot's perspective the body x- 
axis is defined positive forward, the body y-axis 
is positive to the right, and the body z-axis is 
positive downward. Assume that the sensor is 
located at a point P, whose location with respect 
to the reference body axis system is given by the 
position vector r p . The sensor at point P will 
sense the motion (displacement, velocity and/or 
acceleration) at P due to both the rigid aircraft 
and airframe elastic motions. The motion at 
point P on the aircraft in flight is given by the 
linear displacement vector, up, and the angular 
displacement vector, 0p. These displacements at 
point P are obtained by expanding the motion at 
point P in a series of orthogonal free vibration 
modes, in which the first six modes represent 
the rigid body motion: 


oo 


up(rp.t) = £qk(t)5 k(rp) 

(9) 

k=1 


OO 

0p(rp,t) = £qk(t)yk(rp) 
k=1 

(10) 


The orthogonality property of the modes implies 
that the elastic airframe modes produce no net 
displacement of the aircraft center of gravity. 
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The first six degrees of freedom in Eqns. (9) and 
(10) represent the rigid body motions. The 
generalized coordinates qi, q2, and q 3 are the 
rigid body angular motion around the x-axis (roll 
4>, positive right side down), around the y-axis 
(pitch 0, positive nose up), and around the z-axis 
(yaw \|r, positive nose to the right), respectively. 
The generalized coordinates q 4 , qs, and q6 are 
the body linear motions along the x-, y-, and z- 
axis, respectively. The rigid body motion can be 
specified by the linear perturbation velocity of 
the body center of gravity and the angular 
perturbation velocity of the body around the body 
c.g. The linear velocity perturbation of the body 
center of gravity is given by: 


Eqn. (12) 


U C g - (fcg.p x ) RfT 


q2 


Vqay 


(14) 


and 



( OF') 


<qi> 

0p = 

®F 

Ol'Fj 

= Rft 

q2 

U 3 J 


(15) 


where rcg.p is the position vector of point P 
relative to the center of gravity in the body axis 
system, Vft is the perturbation rigid body trim 
Euler angle for yaw, and x indicates a vector 
cross-product. 


Ucg, rigid = 


f : xf^ 

yF 


Therefore the mode shapes for the rigid body 
motion are: 


V* J 


R1 • • • £6] = [ (-rcg.p x)Rft i ] (16) 


and the rigid body angular velocity perturbation 
is given by: 


°>F = [ r ft] 


( dpA 

OF 

VVF ) 


( 12 ) 


M ...76 ] = [Rft 0] (17) 

where 1 represent a 3x3 diagonal unit matrix and 
0 represent a 3x3 zero matrix. 

The total velocity of point P in the body axis 
system is the sum of the trim velocity, Vtrim. 
and perturbation velocities: 


where the subscript F indicates the body axis 
system. The matrix Rft represents the 
transformation of the body axis system to the 
Earth-fixed axis system. This matrix defines the 
aircraft trim attitude with respect to earth axis 
and is given by: 


" 1 0 -sin6FT 

Rft =[ 0 cosqFT sindFTcosGFT 
0 -sinqFT cosqFTcoseFT -I 


(13) 


where 0 ft and 0 ft are the trim Euler angles for 
rigid body pitch and roll. The trim Euler angles 
are obtained from the CAMRAD/JA trim solution 
and are dependent upon the flight conditions. 

The linear and angular displacement perturbation 
vectors of a point P are therefore given by: 


up= Vtrim + X qk 5k 
k=1 


(18) 


The components of Vtrim in the x-, y-, and z- 
directions are obtained from the CAMRAD/JA 
trim solution. The acceleration of point P is the 
sum of the perturbation accelerations and the 
inertial acceleration due to the rotation of the 
trim velocity vector by the body axes angular 
velocity. Therefore in the body axis system: 


up = (of x Vtrim + X qk£k 

k=1 


Eqn. (12) 


(‘Vtrim x ) 


Aqi A 


q2 

Vqay 


+ X qk$k 

k=1 


(19) 


u p = Ucg + 


af^ 

0F 

w 


r cg,p 


Equations (14)-(15), (18), and (19) describe the 
three components of the displacement, velocity, 
and acceleration of point P, respectively, in the 
aircraft body axis system. 
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A separate axis system is defined as being 
attached to the sensor at airframe location P. If 
the orientation of this sensor axis system does 
not coincide with the body axis system the 
displacement vector up, velocity vector Cip, and 
acceleration vector up need to be transformed 
into the sensor axis system by consecutive 
rotations of the vector around the body z-axis, 
the new y-axis, and the resulting x-axis. The 
components of the motion vector in the sensor 
axis system are used to define the elements of 
jth row of the sensor matrices C| (i=0,1,2), which 
in turn define the output of a single-direction 
motion sensor at point P, represented by the j th 
element of the sensor vector y. 

TIME RESPONSE 

The time history response of the linear system 
of Eqn. (8) is given by 16 : 

x s (t) = e A (t-t 0 )xs(to) + 

t 

JeA(t-x)Bv(T)dx (20) 

to 

If the system matrix A has distinct eigenvectors 
the exponential power e^t can be rewritten as: 

eAt = Me( M ' lAM ) 1 M' 1 (21) 

where M is a matrix of linearly independent 
eigenvectors. The columns of matrix M, the so- 
called modal matrix, are made up of the 
eigenvectors with columns corresponding to the 
eigenvalues of matrix A. The matrix 

manipulation M _1 AM results in the 

transformation of A into the diagonal eigenvalue 
matrix a: Thus, Eqn. (20) can be rewritten as: 

x s (t) = Me A (t-t 0 )M-1 x s (to) + 
t 

jMe A (t-T)M-lBv(T)dT (22) 

to 

For a given constant coefficient system, the 
modal matrix, its inverse, and the control 
matrix B are constant coefficient matrices. If 
it is assumed that the control vector v(t) 
(Eqn. (2)) is constant over the considered time 
period from to to t, then Eqn. (22) can be 
rewritten as: 

x s (t) = Me A (t‘to)M ^xs(to) + 


M( J e A(t-T)d T )M‘lBv (23) 

to 

As previously discussed the system control input 
v(t) is defined as: 

v(t)=vp(t)+v ac (t) (2) 

If the control input v(t) is digitized at time step 
intervals of At and the assumption is made that 
the control input is constant over the time step 
At, i.e. v(t)=v(tj) for tj < t < tj+At, then the 
system response xs at time tj+At can be 
calculated from the system state at time tj, 
xs(tj), and the system response to a step input of 
magnitude v(tj) at time tj by using Eqn. (23). 

x s (tj+At) = MeAAtM-1 x s (tj) + 
tj+At 

M JeA(tj-t)dxM-lBv(tj) (24) 

ti 

or in matrix form: 

xs(tj+At) = Me AAt M'1 Xs(tj) + 

A-lM{eAAt.j Bv(tj) (25) 

where I is the diagonal unit matrix. 

The response Xs(tj) in Eqn. (24) represents the 
state of the system at time tj. The state 
velocities *s(ti) are obtained from Eqn. (8). 
Knowledge of the first order state variables xs 
and &s at time tj allows the extraction of the 
displacements x, the velocities Jt, and the 
accelerations X of the aircraft degrees of 
freedom in Eqn. (1). The sensor vector y can then 
be calculated from Eqn. (2) and the active control 
input Vac fro™ Eqn. (3). The input v to the system 
at time step (tj+At) is set equal to ( Eqn. (2)): 

v (tj+At) = v p (tj+At) + Vac(tj) (26) 

NO ISE 

To account for instrumentation signal noise 
and/or for sensor modelling uncertainty within 
the sensor model, noise can be simulated by 
rewriting the sensor vector, y (Eqn. (3)): 

y(tj) = C2X(tj) + Cix (tj) + Cox(tj) + 

Dov(tj) + w(tj) (27) 
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where w(tj) is a vector representing sensor 
noise. For this study the sensor noise is defined 
as: 

w(tj) = rg m , -1.0 <; r £ +1.0 (28) 

where r is a random number and gm defines the 
sensor noise level or magnitude. 

AIRCRAFT MODEL 

An XV-15 size joined-wing tilt-rotor aircraft 
(model 166CL in Ref. 2) in cruise flight was 
considered for this investigation. The aircraft 
has a three-bladed, 25-foot diameter, gimbal- 
mounted, stiff-inplane proprotor. Table 2 
provides a description of the key geometric 
parameters for the rotor and the airframe as 
obtained from Ref. 2 and 1 7. 

The airframe structural characteristics, 
consisting of the modal frequency and mass and 
of the hub motion due to elastic airframe 
deformation for the joined-wing tilt-rotor 
aircraft were reported in Ref. 2 for six airframe 
mode shapes, representing three symmetric and 
three anti-symmetric modes. Table 3 provides 
these data for the six elastic airframe degrees 
of freedom. This mode shape information at the 
hub location is used in the CAMRAD/JA code, 
which provided the mathematical plant model of 
the joined-wing tilt-rotor aircraft. Table 4 
presents the location of the aircraft center of 
gravity and the rotor hubs for the aircraft. 

RESULTS 

The following degrees of freedom were included 
in the present investigation: six rigid-body 
degrees of freedom, six (three symmetric and 
three anti-symmetric) elastic airframe modes, 
two rotor blade bending and one rotor blade 
torsion mode, the hub gimbal, and rotor speed 
degrees of freedom. The plant model, as 
calculated by CAMRAD/JA, consisted of the 
system matrices Ai (i=0,1,2) and the control 
matrix Bo for both symmetric and anti- 
symmetric flight modes. Both the symmetric and 
anti-symmetric plant models consisted of 
eighteen degrees of freedom per flight mode: 
three rigid and three elastic body modes; the 
flap, chord, and torsion modes for each rotor 
blade, resulting in nine rotor modes for the three 
bladed rotor; two hub gimbal and one rotor speed. 
Structural damping was set at 3 percent of 
critical damping for the rotor modes and 2 
percent critical damping for the elastic airframe 
modes. The airframe aerodynamic damping was 
assumed zero. The control vector, v, consisted 


of the collective pitch and the lateral aVid 
longitudinal cyclic pitch. 

STABILITY OF AIRCRAFT WITHOUT FEEDBACK 

This section will discuss the results of the 
stability analysis for the basic aircraft, i.e., the 
joined-wing tilt-rotor aircraft without feedback 
control. The CAMRAD/JA code was used to obtain 
the mathematical plant model of the basic tilt- 
rotor aircraft in cruise as given by Eqn. (1) for a 
forward flight range from 150 to 350 knots at 
standard sea level conditions. For a tilt-rotor 
aircraft configuration in cruise flight the 
CAMRAD/JA analysis divides the aircraft motion 
into uncoupled symmetric and anti-symmetric 
flight modes. Separate system matrices, Aj 
(i=0,1,2) and Bo, are calculated for each flight 
mode. The system matrices for a selected flight 
mode and flight trim condition formed the input 
to the analysis code, which used an eigenvalue 
analysis to evaluate the stability of the aircraft 
without the application of active controls. 

Figure 6 presents the results of the stability 
analysis for the basic joined-wing tilt-rotor 
aircraft in the form of root-locus plots for both 
the symmetric (Figs. 6a and 6b) and anti- 
symmetric (Figs. 6c and 6d) flight modes. The 
natural frequency and damping ratio for the three 
wing modes were calculated and are shown in 
Figs. 7 and 8 for the symmetric and anti- 
symmetric flight modes, respectively. Wing 
modes Is and la are the critical modes for the 
symmetric and anti-symmetric flight modes, 
respectively. The symmetric flight modes shows 
a whirl-flutter velocity of approximately 243 
knots, while for the anti-symmetric flight 
modes the occurrence of whirl-flutter occurs at 
the higher flight speed of approximately 275 
knots. Reference 2 reports a flutter speed of 
approximately 240 and 250 knots for the 
symmetric and anti-symmetric flight modes, 
respectively. The discrepancy between the 
results of the present investigation and of 
Ref. 2 for the anti-symmetric flight mode 
flutter speed is not understood. 

CLOSED LOOP RESULTS - State Variable Domain 

As noted in the previous section the onset of 
whirl-flutter occurred in the symmetric flight 
mode first and therefore the emphasis in the 
present investigation was towards determining 
the feasibility of using active controls to 
increase the flutter velocity for this flight 
condition. The flight velocity of 275 knots was 
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selected as the flight condition for evaluation of 
various active control schemes. 

Accelerometers were considered to be the 
preferred sensor to be used in the feedback 
system. Because of the availability of rotor 
cyclic pitch control on the tilt-rotor aircraft in 
cruise flight it was decided to use the cyclic 
pitch as the active control input v 3 C' 
References 1 0 and 1 1 also used the closed-loop, 
cyclic pitch control to improve the whirl-flutter 
stability of the cantilever-wing tilt-rotor 
models. 

The feasibility of whirl-flutter alleviation 
through active controls was first investigated in 
the state space domain. Only state variable 
accelerations were used in the sensor model. 
The matrices Ci, Co, and Do in Eqn. (3) were zero 
matrices, while the Ca matrix elements had the 
value of "1" or "0", depending upon whether a 
particular state variable acceleration was 
incorporated into the closed-loop feedback 
system. 

Since the first wing mode, q w .is> was the 
critical mode for the onset of whirl-flutter, the 
acceleration for this state variable was selected 
for the sensor output and the longitudinal cyclic 
pitch, 5 S , and lateral cyclic pitch, 5 C , were 
selected for possible active control inputs. 

Figure 9 presents the effect of various 
longitudinal cyclic pitch feedback gains on the 
stability of the elastic wing mode Is at 275 
knots. A positive feedback gain for q w ,is/8s has 
the desired stabilizing effect on the wing mode 
Is, has little effect on wing mode 3s, but has a 
destabilizing effect on wing mode 2s. 
Figure 10 shows the effect of lateral cyclic 
pitch feedback on the frequency and damping 
ratio of wing modes Is, 2s, and 3s at 275 knots. 
A negative feedback gain for 3w,is/S c has a 
stabilizing effect on wing mode Is, while having 
a destabilizing effect on wing modes 2s and 3s. 
Figures 9b and 10b show that all three wing 
modes will be stable for 1.6<“q w ,is/S s <4.0 g/rad 
or q w .is/8 c <-2.5 g/rad. 

It should be noted here that feeding back the 
wing mode Is acceleration as described above 
also influences the frequency and damping of 
other body and rotor degrees of freedom, i.e., 
their eigenvalue location in the root locus plot 
changes. One rotor mode, a rotor lead/lag 
bending mode, was significantly influenced by 
the feedback of wing mode Is acceleration, 
qw.is. to either the longitudinal or lateral cyclic 
pitch. This sensitivity is shown in Fig. 11. 


From comparison of Figs. 9b and 11b it is se^n 
that at the feedback gain values of 
1 ,6< "qw,i s /8 s <4.0 g/rad where the wing mode Is 
is damped this critical rotor mode has become 
unstable. Comparison of Figs. 10b and 11b 
shows that this rotor mode has also become 
unstable at feedback gain values for 
qw.is/5 c <-2.5 9 /ra d, which is the desired gain 
level to stabilize wing mode Is. Figure 11b 
shows that the magnitude of negative damping 
ratio for the critical rotor mode is slightly 
lower for the required feedback gain levels for 
qw.is/5 s as compared to the required feedback 
gain levels for q w .is/$c t0 stabilize the wing 
mode Is. Therefore the longitudinal cyclic pitch, 
5 S , was selected as the active control input. 

An obvious choice to stabilize the system is the 
additional feedback of the critical lead/lag mode 
to cyclic pitch. A sensitivity analysis showed, 
however, that this critical rotor mode could be 
influenced by feedback of wing modes 2s and 3s 
to either longitudinal or lateral cyclic pitch. 
Again the longitudinal cyclic pitch, 5 S , was 
selected at the active control input for feedback 
of wing modes 2s and 3s. 

As this study was a preliminary investigation 
into the use of active controls for whirl-flutter 
alleviation, no effort was made to optimize the 
feedback control system. A number of 
combinations of cjw,is/8s> fW,2s/Ss« <ind €|w,3s/8s 
gain values were evaluated using the closed-loop 
analysis code until a stable system (including 
rotor stability) was obtained for the desired 
condition of 275 knots flight velocity. 
Satisfactorily stability was obtained for a 
closed-loop system using a qw.is^s gain of 2.0 
g/rad and a cjw,3s/S s gain of -1.25 g/rad. The 
effect of flight velocity on the stability of this 
closed-loop symmetric flight mode system is 
shown in Fig. 12. The natural frequencies and 
damping ratios for the three wing modes for both 
the basic aircraft and the closed-loop system 
are shown. Little effect of using this type of 
active control on the natural frequency of these 
wing modes is observed (Fig. 12a). However, the 
damping ratio of all three wing modes is changed 
substantially (Fig. 12b). Wing mode Is shows 
more damping, while wing modes 2s and 3s show 
reduced damping from the configuration with no 
feedback control. Figure 12 shows that even 
non-optimal feedback in the state domain of 
wing chord and wing torsional accelerations to 
longitudinal cyclic pitch can increase the tilt- 
rotor whirl-flutter speed from 243 to above 275 
knots, a 13.2% improvement in the flutter 
velocity. 
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CLOSED LOOP RESULTS - Physical Domain 

The previous section showed that feedback of 
state variable accelerations to cyclic pitch can 
be used to delay the onset of whirl-flutter on the 
joined-wing tilt-rotor aircraft. However, to be 
able to feed back pure state variables requires 
the identification of the individual system 
modes which would potentially require a large 
number of sensors mounted on the aircraft. 
Flutter alleviation studies on fixed-wing 
aircraft have shown that it is possible to 
increase the aircraft flutter velocity by feeding 
back the output of a limited number of 
accelerometers mounted at appropriate locations 
near the wing tip to an active control input, such 
as a wing control surface 18 - 19 . The use of a 
limited number of accelerometers mounted at 
appropriate locations near the wing tip, whose 
signal output is then fed back to a blade pitch 
control was studied in Ref. 11 for a cantilever- 
wing tilt-rotor aircraft configuration to 
increase the whirl flutter speed from 285 to 
above 300 knots. Use of a limited number of 
wing mounted accelerometers was also 
investigated here for the joined-wing tilt-rotor 
aircraft. Again the rotor longitudinal cyclic 
pitch was selected for the active control input. 

The results of the state variable domain analysis 
discussed in the previous section showed that 
the feedback of the accelerations of the wing 
modes Is and 3s degrees of freedom to 
longitudinal cyclic pitch could be used to 
stabilize the joined-wing tilt-rotor aircraft at 
275 knots flight speed. It was therefore decided 
to use wing-mounted, single-direction 
accelerometers which would measure these two 
aircraft degrees of freedom as the sensors in the 
feedback system. Only the symmetrical flight 
mode was analyzed in the "physical" domain. In 
application, wing-mounted sensors would also 
measure the acceleration due to the anti- 
symmetric wing elastic deformations and anti- 
symmetric rigid body motion. However, by 
placing sensors in the same location on opposite 
wings and by summing the corresponding sensor 
output signals the contribution of these anti- 
symmetric flight modes could be eliminated. 

As mentioned earlier, Ref. 2 only provided the 
elastic airframe mode shape information for the 
rotor hub location. Therefore hub-mounted 
accelerometers, which provide hub accelerations 
in the body axis system were considered here. 
Reference 1 1 showed that the location of the 
wing-mounted sensors can be an important 
parameter in the closed loop active control 
system to increase the whirl flutter velocity on 
a tilt-rotor aircraft. Therefore the hub location 


is not necessarily the best location for ftie 
closed-loop system sensors. 

The airframe modal information for the hub 
location together with the CAMRAD/JA 
calculated trim information (the cruise velocity 
vector Vtrim and the trim Euler angles 4>ft and 
6ft, Eqn. 11) were input to the sensor model 
utility (Fig. 4). This utility then determined the 
contribution of each aircraft degree of freedom 
to the hub accelerations in body x-, y-, and z- 
directions. The assumption was made that the 
contribution of the symmetric rigid body motion 
to the acceleration at the hub (Eqn. (17)) was 
eliminated by means of a high pass filter. 
Mathematically this equates to setting the 
corresponding Cj-matrix elements to zero. As 
mentioned previously, no attempt was made to 
model the dynamics of the swashplate actuators 
in the mathematical description of the feedback 
control system. The 275-knot flight trim 
condition was again selected for the evaluation 
of the various feedback control schemes. 

Table 3 shows that for wing mode Is the linear 
displacement in vertical direction is an order of 
magnitude higher than the linear displacements 
in both the chord-wise and span-wise directions. 
A single-direction, hub-mounted accelerometer 
with its sensing direction along the body z-axis 
was therefore selected for the measurement of 
the wing mode Is degree of freedom. As was 
expected from the results of the state domain 
analysis, negative feedback of the hub vertical 
acceleration, Uh,v> to the longitudinal cyclic 
pitch, S s , also adversely affected the stability of 
the critical rotor lead/lag mode (at 2.3/rev) as 
seen in Fig. 13a. A sensitivity analysis showed 
that this adverse effect could be reduced by 
eliminating the contributions of wing modes 2s 
and 3s to the hub vertical acceleration 
measurement as can be seen from comparison of 
Figs. 13a and 13b. Physically this could be 
accomplished be means of a notch filter on the 
accelerometer output signal. Mathematically the 
effect of the notch filter was simulated by 
setting the C 2 -elements' corresponding to the 
mode 2s and 3s contributions to zero. 
Eliminating the wing mode 2s and 3s 
contributions also reduces the effect of the 
Uh.v/8s feedback on the damping ratios of these 
two wing modes as seen from comparison of 
Figs. 13a and 13b. 

A sensitivity analysis showed that feeding back 
the hub acceleration in wing chord-wise 
direction, iih.c. and in wing span-wise direction, 
Uh.s. to longitudinal cyclic pitch, 5 S , had a 
stabilizing effect on the critical rotor mode as 
shown in Figs. 14a and 14b, respectively. A 
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slight beneficial effect of such feedback on the 
Is wing mode damping ratio is also observed in 
Figs. 14a and 14b, although the mode remains 
unstable. 

A feedback system was selected consisting of 
two hub accelerometers measuring vertical and 
span-wise hub accelerations, whose outputs 
were fed back to the longitudinal cyclic pitch. 
The hub vertical acceleration measurement was 
assumed to be conditioned by means of a notch 
filter so as to eliminate the contribution of wing 
modes 2s and 3s from this sensor output signal. 
As mentioned before a high pass filter was also 
used to eliminate the contributions of the rigid 
body modes to the measured hub accelerations. 
This feedback system was analyzed in the 
following manner. The value of the feedback 
control gain GhV$s was fixed while the feedback 
gain Uh.s/Ss was varied systematically within the 
closed-loop analysis code. After each Uh,s/5 s 
change the code performed an eigenvalue 
analysis to determine the stability of the 
closed-loop system. This process was then 

repeated for different Gh,v/6 S values. 

Feedback of hub vertical acceleration Cih.v 
(gain=0.0048g/rad, with a notch filter on the 
sensor output signal), and hub span-wise 
acceleration Gh.s. (gain=-0.0035g/rad) to the 

longitudinal cyclic pitch produced a stable 

closed-loop system. This feedback system was 
analyzed over the tilt-rotor aircraft velocity 
range from 150 to 350 knots. Figure 15 
compares the stability of the three wing modes 
for the basic joined-wing tilt-rotor aircraft (no 
feedback) and for the closed-loop system. Little 
effect on wing natural frequencies is found 
(Fig. 15a). The effect of flight speed on the 
damping for the closed-loop system is shown in 
Fig. 15b. Comparison of Fig. 15b with the 
results obtained for the state variable domain 
analysis (Fig. 12b) shows similar results for the 
damping ratio for the critical wing mode Is. 
Lower damping ratio values are observed for 
wing mode 2s for the closed-loop "physical" 
domain analysis (Fig. 15b) when compared to the 
state variable domain results (Fig. 12b). A large 
decrease in the wing mode 3s damping ratio was 
seen in Fig. 12b for the closed-loop system as 
compared to the no feedback case. In contrast, 
Fig. 15b shows very little influence of the 
"physical" domain closed-loop feedback system 
on the wing mode 3s damping ratio. Figure 15 
illustrates that again a 13.2% increase in the 
whirl-flutter velocity from 243 to 275 knots 
was obtained. 

Figure 16 shows the root locus plot for this 
closed-loop system for the speed range from 150 


to 350 knots. Comparison of the root locus plSts 
of Fig. 6 (basic aircraft without feedback 
control) with Fig. 16 (closed-loop system) also 
shows the improved stability for the closed-loop 
system, but at the same time illustrates the 
effect of the feedback system on the damping of 
the various rotor and aircraft degrees of 
freedom. The rotor lead/lag mode at 2.3/rev 
shows an decrease in damping due to the hub 
acceleration feedback system and is neutrally 
stable at approximately 278 knots. Little effect 
on the other rotor modes and on the rigid body 
modes is observed. 

TIME RESPONSE CALCULATIONS: 

To evaluate the magnitude of the required active 
control input, the time response of the aircraft 
at a 270-knot cruise condition was calculated. 
Time-history response calculations were made 
for both the basic aircraft (no feedback control) 
and for the closed-loop system. A longitudinal 
cyclic pitch control input, v p (Eqn. (2)), 
representing a square doublet with a 1.0 deg 
magnitude and a period of 1.0 sec (Fig. 17a) was 
used to excite the system at this 270-knot trim 
condition. The calculated sensor response for 
the two accelerometers at the rotor hub is 
presented in Figs. 17b and 17c for the basic 
aircraft. The response of the higher aircraft 
modes results in the higher frequency response 
in the accelerometer output (Figs. 17b and 17c) 
during the first 1-1.5 sec (approximate 12 rotor 
revolutions) of the motion. A rapid divergence of 
the vertical acceleration magnitude is observed 
(Fig. 17c). The acceleration in the vertical 
direction is approximately a factor 5 higher than 
the measured span-wise acceleration (Fig. 17c 
vs. Fig. 17b). 

The time response of the closed-loop system to 
the same square doublet is shown in Fig. 18 for 
the first 2 seconds and in Fig. 19 for the first 
15 seconds, which equates to approximately 115 
rotor revolutions. The closed-loop system 
(Fig. 18b) shows slightly higher span-wise 
accelerations during the first 2 seconds of the 
response as compared to the basic aircraft time 
response (Fig. 17b). The vertical accelerations 
for the closed-loop system (Fig. 18c) are 50% 
lower than those observed for the basic aircraft 
(Fig. 17c) during this 2 second period. 
Figure 19 shows that after 15 seconds (115 
rotor revolutions) the aircraft has not yet 
reached a steady-state condition. The measured 
accelerations in vertical and span-wise 
direction after 15 seconds are approximately 
one-third (0.5g's and 0.2g's, respectively) the 
maximum value experienced due to the doublet 
excitation (1.9g’s and 0.7g’s, respectively). 
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After 10-12 seconds the hub vertical 
acceleration measurement shows a single 
frequency content of 3.6 Hz. This corresponds to 
the first wing mode (mode Is). The time 
response signal for the span-wise accelerometer 
is presented in Figs. 18b and 19b and shows a 
multiple frequency content throughout the 15 
seconds of the calculated aircraft time response. 
Again a higher frequency response is observed in 
the accelerometer outputs (Figs. 18b and 18c), 
which is similar in character to the sensor 
output of the basic aircraft (Figs. 17b and 17c) 
during the first 1.5 seconds. It is believed that 
the higher frequency content during the first few 
seconds is a basic aircraft behavior and is not 
the result of the feedback system. However, the 
influence of these higher frequencies on the 
measured signal persist for approximately 10 
seconds (Figs. 19b and 19c), and it is thought 
that the strong coupling of the rotor lead/lag 
mode at 2.3/rev with the first wing mode (Is) 
within the closed-loop system is the main cause 
of the longer time period for the higher 
frequency influence to damp out. 

The longitudinal cyclic pitch active control input 
is shown in Figs. 18a and 19a. A maximum 
active control input of 0.009 deg is observed. 
The 0.009 deg cyclic pitch input equates to a 
rotor active control force of approximately 9 lbs 
being applied at the rotor hub to damp the wing 
motion and to delay the tilt-rotor whirl-flutter 
instability. This level of active control input for 
the longitudinal cyclic pitch is of the same 
magnitude as that predicted by the author for a 
similar flutter alleviation study for a 

cantilever-wing tilt-rotor aircraft 
configuration 11 and is a major improvement 
over the levels of cyclic pitch control predicted 
by Nasu 7 , who showed a required cyclic pitch 
control angle of 0.5 deg. The higher frequency 
content can also be observed in the active 
longitudinal cyclic control input in Fig. 18a and 
in Fig. 19a for the first 10 seconds of the 
motion. After 10-12 seconds the active control 
input shows the single frequency content of 
3.6 Hz, corresponding to the critical wing 
mode Is. 

For the actual vehicle it is unlikely that 
actuators can be accurately controlled for 0.009 
deg of pitch input. Control system dynamics 
(actuators, swashplate, hydraulic system,...) will 
have to be included in the simulation to achieve 
more realistic results. However, these results 
provide two important findings. First, the level 
of hub force generation necessary to stabilize 
the tilt-rotor is very small and can be easily 
withstood without adverse dynamic loads or 
degrading aircraft handling quality 


characteristics. Second, the frequerfcy 
requirements for the control actuators are 
within the bandwidth of current electro- 
hydraulic actuation flight hardware. 

The closed-loop system response to a square 
doublet was also calculated for two cases where 
the sensor signals were contaminated by noise, 
representing signal instrumentation noise and/or 
uncertainty regarding the accuracy of the system 
modeling. Based upon the time response 
calculations shown in Fig. 19, noise gain levels 
for g m (Eqn. (28)) of 0.1 and 0.5 were selected. 
This represents a noise level of 0.1 g and 0.5g, 
respectively, for the two hub-mounted 
accelerometers. The corresponding time 
responses for the closed-loop system are shown 
in Figs. 20 and 21 , respectively. Comparison of 
Fig. 19 (no noise) with Fig. 20 (0.1 g noise) 
shows essentially the same time response 
behavior with the noise contamination. 
Figure 21a shows that the 0.5g sensor noise 
contamination level causes an increase in the 
required level of active longitudinal cyclic pitch 
control from 0.009 deg (no noise) to 0.010 deg 
for the first few seconds after disturbing the 
system. After 10 seconds an increase in the 
required magnitude of the active control level 
from 0.0025-0.0030 deg for the no-noise case 
(Fig. 19a) to approximately 0.005 deg for the 
0.5g noise case is observed (Fig. 21a). The 
periodic character of the active control is not 
readily observable in Fig. 21a. A decrease in the 
magnitude of the hub horizontal and vertical 
accelerations over time is seen (Figs. 21b and 
21c). The observed response decay rate for the 
0.5g noise case (Fig. 21 ) is smaller than the 
decay rates observed for the no-noise (Fig. 19) 
and 0.1 g noise (Fig. 20) cases. However, the 
closed-loop system with 0.5g noise 
contamination on the sensors, representing 25% 
of the maximum observed sensor response level 
of 2.0g vertical hub acceleration, still shows a 
stable system indicating a robust control 
system. 

CONCLUSIONS 

The present investigation showed that it is 
possible to delay the onset of whirl-flutter on a 
XV-15 size joined-wing tilt-rotor aircraft 
through the use of active controls. The 
CAMRAD/JA code was used to obtain a set of 
linear differential equations, which describe the 
motion of the tilt-rotor aircraft in cruise. The 
airframe vibration mode shapes, modal 
frequencies, and modal masses for a 
representative joined-wing tilt-rotor aircraft 
were used. The CAMRAD/JA output, consisting of 
the open-loop system matrices, and a sensor 
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model utility output were input to a separate 
program, which performed the closed-loop, 
active control calculations. An eigenvalue 
analysis was performed to determine the flutter 
stability of both open- and closed-loop systems. 

• Sensor models based upon the feedback of 

pure state variables increased the whirl- 
flutter velocity from 243 to over 275 knots 
by feeding back the first and third wing 
elastic bending degrees of freedom to the 
longitudinal cyclic pitch. 

• Whirl-flutter could also be delayed from 243 

to above 275 knots by feeding back hub span- 
wise and vertical accelerations to the 

longitudinal cyclic pitch. 

• Time response calculations at a 270-knot 
cruise condition showed an active cyclic pitch 
control level of 0.009 deg, which equates to a 
very acceptable 9 pound active control force 
applied at the rotor hub. 

• Contamination of the sensor output signal 

with noise, whose magnitude was equal to 

25% of the maximum measured amplitude of 
the sensor output did not adversely effect the 
closed-loop system stability, indicating a 
robust control system. 
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Appendix A; Comprehensive Rotorcraft Analysis 

CAMRAD/JA (for Comprehensive Analytical Model 
of Rotorcraft Aerodynamics and Dynamics, 
Johnson Aeronautics version) is an analysis 
designed to calculate rotor performance, 
aerodynamic and structural loads, aircraft 
vibration, gust response, flight dynamics, 
handling qualities, and aeroelastic stability. The 
analysis development is discussed in detail in 
References 12-15. 

The rotor aerodynamic model is based on lifting- 
line theory, and uses two-dimensional airfoil 
characteristics, and a vortex wake. For the 


2-12 



aeroelastic stability analysis it is generally 
sufficient to use the uniform inflow model of 
CAMRAD/JA, i.e. assuming a linear variation of 
inflow over the rotor disk. The rotor structural 
model is based on engineering beam theory for 
rotating wings with large pitch and pretwist. 
Both rigid and elastic blade motion are included. 
A modal method is used to calculate, in vacuum, 
the blade bending modes for the rotor. Blade 
torsion modes are not coupled with the bending 
modes, but are calculated independently. The 
total elastic rotor response is a linear 
combination of the modal solution in which the 
bending and torsion motion of the rotor are 
coupled. 

An orthogonal mode representation of the body 
elastic motion is used in CAMRAD/JA. The 
airframe structural vibration modes need to be 
obtained from a separate structural analysis (in 
the present analysis the MSC-PAL finite element 
structural analysis code was used). The 
linearized equations of motion for the rigid body 
degrees of freedom are used in the calculation of 
the aircraft vibration and transition motion. For 
this study the aerodynamic forces on the wing- 
body, horizontal tail and vertical tail are 

modelled, including control surfaces. Static 
aircraft aerodynamic characteristics were 
obtained from a wind tunnel test. The 
generalized aerodynamic damping and control 
forces on the airframe elastic modes can be 
included in CAMRAD/JA and these terms are 
normally estimated for the frequency of the 
principal excitation of the mode. In the present 


analysis these dampings were set to zero, 'ftie 
airframe aerodynamic representation is 
basically a quasi-static model. 

The aircraft motion consists of the six rigid- 
body degrees of freedom and the elastic free 
vibration body modes. A body-axis coordinate 
frame with origin at the aircraft center of 
gravity is used for the description of the motion. 
For the elastic motion of the aircraft in flight, 
the displacement and rotation at an arbitrary 
point are expanded in a series of the orthogonal 
free vibration modes. The first six degrees of 
freedom for the elastic aircraft motion at this 
arbitrary point are the rigid body motions, while 
the remaining body degrees of freedom represent 
the elastic modes of the aircraft. (See the 
sensor model derivation in the main text). 

The rotor equations of motion require the six 
components of the hub linear and angular motion 
in the shaft axis system due to both the rigid and 
elastic aircraft body motion. This information 
was obtained from the results of the MSC-PAL 
analysis as reported in Ref. 2 (see table 3). 

CAMRAD/JA first performs the trim analysis, in 
which the equations of motion are solved for the 
case of a steady state flight condition. 
Aeroelastic stability is then calculated from the 
trim solution by constructing the system 
matrices that describe the linear differential 
equations of motion and performing an 
eigenvalue analysis. 
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Table 1: Comparison of tilt-rotor model of Refs. 10 and 11 with present investigation. 


< 


Nasu (Ref. 10) 

van Aken (Ref. 11) 

Present investigation 

Cantilever-wing model 

Cantilever-wing aircraft 

Joined-wing aircraft 

Semi-span wing model: 

Free flight: 

Free flight: 

no rigid body motion 

rigid body motion 

rigid body motion 

symmetric flight mode 

symmetric and anti-symmetric 

symmetric and anti-symmetric 


flight modes 

flight modes 

Elastic body modes: 

Elastic body modes: 

Elastic body modes: 

symmetric modes only 

symmetric and anti-symmetric 

symmetric and anti-symmetric 

wing chordwise bending 

wing chordwise bending 

wing mode 1 

wing beamwise bending 

wing beamwise bending 

wing mode 2 

wing torsion 

wing torsion 

wing mode 3 


pylon yaw 


Rotor: 

Rotor: 

Rotor: 

3 bladed hingeless rotor 

3 bladed hingeless rotor 

3 bladed hingeless rotor 

blade chordwise bending 

blade chordwise bending 

blade chordwise bending 

blade beamwise bending 

blade beamwise bending 

blade beamwise bending 


blade torsion 

blade torsion 


rotor speed 

rotor speed 


gimbal motion 

gimbal motion 

Feedback system: 

Feedback system: 

Feedback system: 

using state variables 

using state variables 

using state variables 


using wing mounted 

using hub mounted 


accelerometers 

accelerometers 


Table 2: Specifications of joined-wing tilt-rotor aircraft model 


Rotor: 


number of blades 

3 

blade radius 

12.5 ft 

solidity, 0 

0.0890 

airfoil section 

64-series 

rotational speed 

458 rpm 

shaft cant angle 

-1 deg. 

mast height 

4.667 ft 

Wing: 


semispan 

16.1 ft 

forward wing: 


chord 

3.5 ft 

sweep (forward) 

6.5 deg 

airfoil thickness 

12% 

aft wing: 


chord 

1.7 ft 

sweep (forward) 

15 deg 

airfoil thickness 

12% 

airfoil section 

64A212 

Aircraft weight 

13,000 lbs 
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Table 3: Wing mode shape information at right hub 


< 


symmetric modes 

mode Is 

mode 2s 

mode 3s 

frequency, Hz 

3.65 

5.56 

7.30 

modal mass, slug/12 
linear displacements: 

1.833 

0.832 

0.375 

x, in 

-0.0735 

0.1265 

0.0341 

y, in 

0.0771 

-0.1318 

-0.0340 

z, in 

-0.4530 

-0.1439 

-0.0655 

angular displacements: 




x, rad. 

-0.00313 

-0.00001 

0.0015 

y, rad. 

0.00515 

0.00400 

0.00178 

z, rad. 

0.0011 

-0.00199 

-0.00049 


anti-symmetric modes 

mode la 

mode 2a 

mode 3a 

frequency, Hz 

4.75 

7.42 

8.13 

modal mass, slug/12 

0.639 

0.393 

1.295 

linear displacements: 




x, in 

-0.0025 

0.0324 

0.1576 

y, in 

-0.0089 

-0.0104 

-0.1706 

z, in 

0.2578 

-0.0841 

0.0339 

angular displacements: 




x, rad. 

0.001 1 

0.0017 

0.0010 

y, rad. 

-0.0053 

0.00253 

-0.0006 

z, rad. 

0.0001 

-0.00048 

-0.00235 


Table 4: Location of c.g. and rotor hubs on joined-wing tilt-rotor aircraft. 



Body 

axis system coordinates 


X 

(inches) 

y 

(inches) 

z 

(inches) 

Component location: 
center of gravity 

0 . 

0 . 

0 . 

right hub 


192.99 

-26.44 

left hub 

'Hr ft 

-192.99 

-26.44 
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Fig. 3: Effect of nacelle c.g. location on the flutter speed of the 
joined- wing tilt-rotor aircraft (baseline cantilever-wing 
aircraft flutter speed shown for comparison). 

(source: Ref. 2) 
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Fig. 4: Methodology of whirl-flutter alleviation study 



Fig. 5: Body axis system 
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(b) exploded view 
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Fig. 6: Root locus plots for basic joined-wing aircraft without 
active controls. 

(a) symmetric flight mode 

(b) symmetric flight mode exploded view 

(c) anti-symmetric mode 

(d) anti-symmetric mode exploded view 
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Fig. 7: Natural frequency and damping ratio of wing modes for 
basic aircraft without active controls, symmetric flight 
condition; 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 
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Fig. 8: Natural frequency and damping ratio of wing modes for 
basic aircraft without active controls, anti-symmetric flight 
condition; 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 
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Feedback gain q w ls /8 s , g/rad 

Fig. 9: Natural frequency and damping ratio for wing modes 
for closed loop system (symmetric flight mode) when 
feeding back wing mode state variable acceleration, 
q w ls , to longitudinal cyclic pitch, 8 s , at 275 knots 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 



Fig. 10: Natural frequency and damping ratio for wing modes 
for closed loop system (symmetric flight mode) when 
feeding back wing mode state variable acceleration, 
q w Js , to lateral cyclic pitch, 5 c , at 275 knots 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 
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Damping ratio, £ Natural frequency, Hz 



Feedback gain, g/rad 

Fig. 11: Natural frequency and damping ratio for critical rotor 
mode for closed loop system (symmetric flight mode) 
when feeding back wing mode, q w j acceleration to 

longitudinal (8p or lateral (8^ cyclic pitch at 275 knots. 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 



150 200 250 300 350 

Flight Speed V , knots 

Fig. 12: Natural frequency and damping ratio for wing modes 
for basic aircraft and for closed loop system using 
sensors measuring state variable accelerations (symmetric 
flight mode): *q' w ^ = 2.0 g/rad, q* w 3s /S s = -1.25 g/rad 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 
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Fig. 13: Damping ratio for wing and critical rotor modes for 

closed loop system (symmetric flight mode) when feeding 
back vertical hub acceleration to longitudinal cyclic 
pitch. 

(a) without a notch filter on vertical accelerometer output 

(b) with a notch filter on vertical accelerometer output 
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Fig. 14: Damping ratio for wing and critical rotor modes for 
closed loop system (symmetric flight mode) with 
feedback to longitudinal cyclic pitch. 

(a) feeding back chord-wise hub acceleration, *1^ c 

(b) feeding back span- wise hub acceleration,*^ ? 
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Natural frequency, Hz 




Rg. 15: Natural frequency and damping ratio for wing modes 
for basic aircraft and for closed loop system using 
sensor at hub measuring "physical" accelerations: 
u /8 = 0.0048 g/rad, u /8 = -0.0035 g/rad 

• h,» » 

(a) Natural frequency, Hz 

(b) Damping ratio, £ 
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Fig. 16: Root locus plots for closed -loop system using feedback of 
vertical and spanwise accelerations at the hub to 
longitudinal cyclic pitch. 

(a) symmetric mode 

0>) symmetric mode exploded view 
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Spanwise acceleration at hub, g Control input 8 s , deg. 


Elapsed time, sec 

Fig. 17: Basic aircraft time response to a square doublet 8 S 
longitudinal cyclic pitch input; 

(a) Active control input 6 s , deg 

(b) Span-wise acceleration ^ g 

(c) Vertical acceleration ^ y . g 
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Elapsed time, sec 

Fig. 18: Closed-loop time response to a square doublet 8 


(a) Active control input 8 s> deg 

(b) Span-wise acceleration g 

(c) Vertical acceleration*^ g 
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Fig. 19: Closed-loop time response to a square doublet 5 s 
input. 

(a) Active control input 5 g , deg 

(b) Span -wise acceleration $l g 

(c) Vertical acceleration g 
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Fig. 20: Closed-loop time response to a square doublet 8 g 


input with a O.lg sensor noise contamination. 

(a) Active control input 8 gt deg 

(b) Span- wise acceleration u^ g , g 

2-25 (c) Vertical acceleration u^ vt g 





Elapsed time, sec 

Fig. 21: Closed-loop time response to a square doublet 8 

input with a 05g sensor noise contamination. 

(a) Active control input 8^ deg 

(b) Span-wise acceleration uj^.g 

(c) Vertical acceleration y , g 







